Crossover from Fermi-Pasta-Ulam to normal diffusive behaviour in heat conduction 
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We study heat conduction in one, two and three dimensional anharmonic lattices connected to 
stochastic Langevin heat baths. The inter-atomic potential of the lattices is double-well type, i.e., 
Vdw(k) = k2X^ /2 + k4,x'^ /A with fe < and kt, > 0. We observe two different temperature regimes 
of transport: a high-temperature regime where asymptotic length dependence of nonequilibrium 
steady state heat current is similar to the well-known Fermi-Pasta-Ulam lattices with an inter- 
atomic potential, Vp-pvix) = k2X^ /2 -\- kix'^/i with k2,ki > 0. A low temperature regime where 
heat conduction is diffusive normal satisfying Fourier's law. We present our simulation results at 
different temperature regimes in all dimensions. 

PACS numbers: 44.10.-l-i, 05.40.-a, 05.60.-k, 05.70.Ln 



Search for a microscopic lattice Hamiltonian showing 
diffusive normal heat conduction satisfying Fourier's law 
has been an open long-standing problem (Ij, i^] . It is now 
widely accepted that heat conduction in one-dimensional 
(ID) momentum-conserving lattice models is anomalous 
and does not obey Fourier's law, J{x) = —kWT{x), relat- 
ing the local heat current density J(x) at a point x to the 
local temperature T(x) gradient (k defines the thermal 
conductivity which is expected to be an intrinsic property 
of the material) [3| . An exception is a chain of coupled 
rotators exhibiting normal transport properties even in 
the absence of an on-site potential (i.e., momentum con- 
serving lattice) [J, |5| . This model shows a transition from 
infinite to normal thermal conductivity with the increase 
of the local temperature. The typical Hamiltonian of 
a homogeneous ID momentum conserving lattice with 
inter-site potential V{x) is 
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where xi is a scalar displacement at lattice site I. We 
use fixed boundary, i.e., xq = xn+i = 0. The inter- 
site potential for a chain of coupled rotators is given 
by V{x) = 1 — cos(a;) where x is relative rotation be- 
tween two neighboring sites of the chain. A variant of 
the coupled rotator model, namely a double-well poten- 
tial Vuw{x) = k2X^/2 + kix'^/A with ^2 < and k^ > 
has been also studied in Ref.jjj. It is claimed that heat 
conduction in this model shows normal behaviour at low 
temperature. In the first part of this paper we study the 
double-well (DW) potential model in ID rigorously at 
low and high temperatures using simulation with classic 
Langevin model of heat baths. We observe two different 
temperature regimes: a high-temperature regime where 
asymptotic length (N) dependence of thermal conductiv- 
ity is K ~ iV^/"^ in ID. This asymptotic length dependence 
is well-known for ID Fermi-Pasta-Ulam (FPU) lattices 
with an inter-site potential Vfpu(2;) = k2X^ jl + ki^x'^ jA, 



FIG. 1: A schematic of the double- well inter-particle potential 
showing height of the potential barrier k^T\ — k^jAk^. 



where A:2 , ^4 > [6| . One interesting feature of the DW 
model is that wc find the asymptotic length dependence 
of K at high temperature for relatively much smaller 
length (almost two order of magnitude) compared to the 
length scale where it has been shown for the FPU chains 
iTI]. This occurs because of the absence of low energy 
long wavelength heat carrying modes in the DW model. 
We could not unambiguously confirm normal heat con- 
duction at low temperature in this model. We present 
detailed simulation results in this regime and give sev- 
eral arguments supporting our results. 

Next we consider lattices with DW type inter-atomic 
potential in higher dimensions. From our results for ID 
we expect the asymptotic length dependence for higher 
dimensions at high temperature to be similar to that of 
the FPU lattices. It has been predicted that k ^ log(Af) 
for a two dimensional (2D) FPU lattice where A'' is length 
along the applied temperature difference [6, 8] . Although 
available numerical simulations for 2D FPU system show 
K - A^°-22 [i-[i3|. For a three dimensional (3D) FPU 
lattice, based on kinetic theory and Boltzmann equation 
approaches, the expectation is that the thermal conduc- 
tivity should be finite (and Fourier's law is valid) at high 
temperature where Umklapp processes dominate. Recent 
works 



111 Il2| claimed to see evidence of normal heat 



conduction for the 3D FPU lattices even in the absence 
of on-site pinning potential. We find in our simulation 
K ^ log(A^) for the 2D DW model at high temperature. 
Thus we give first confirmation of the predicted logarith- 
mic behaviour for the 2D FPU lattices from our sim- 



ulation. The high-temperature heat conduction in the 
3D DW model does show signature of normal diffusive 
behaviour. We also simulate low-temperature thermal 
transport for the DW model in higher dimensions, and 
we find normal diffusive heat conduction showing a finite 
K for 2D and 3D lattices. 

The DW model was studied in Ref.[4] using Nose- 
Hoover model of thermostat which does not ensure cor- 



rect energy distribution for the system [13| . We here use 
stochastic Langevin heat baths and ensure local thermal 
equilibration in our simulation. The DW model experi- 
ences a structural phase transition with the reduction of 
its temperature, and it shows up in different asymptotic 
length dependence of heat current at low and high tem- 
peratures. We can have a rough estimate of the crossover 
temperature from the following arguments. The height 
of potential barrier between the stable minima and the 
unstable maxima is given by AC/ = fc|/4fc4 (see FigllJ. 
Comparing AC/ with the mean temperature of the baths 
we get a crossover temperature Ti scale for the ID DW 
model, i.e., fceTi = kj/Ak^ with t = (Tl + Tn)/2. In 
our simulation we connect the two end sites of the chain 
described by Eq. Q along with the DW form of the inter- 
site potential to two Langevin heat baths which are kept 
at temperatures Tl and Tr. The bath / introduces a 
fluctuating zero-mean Gaussian white noise rji and a dis- 
sipative term 7; in the equation of motion of the end 
particle of the chain. The terms 7; , rji are related by 
fluctuation-dissipation theorem as given in EqlS) The 
coupling strength to the bath is controlled by the dissi- 
pation constant 7; . The equations of motion for the sites 
on the wire are: 

xi = ~k2{2xi ~ xi^i ~ xi+i) + ki{xi+i - xi)^ 
- k4{xi - xi^if + {--fi^xi + rii,)Si^i 
+ {-IRXN +'r]R)^i.N, (2) 

where the noise-noise correlation is given by 

imithmit')) = 2^ikBTiSlmS{t-t') . (3) 

We set the Boltzmann constant ks — 1- We simulate 
the above set of equations of motion using a velocity- 
Verlet algorithm. We measure the heat current and local 
temperature in the nonequilibrium steady state (after 10^ 
to 10^ time steps of transient dynamics) using following 
definitions. We define time-averaged local heat current 
at site I as 

1 /■* 

f = {ji{t))t=\imt^^- / ji{T)dT (4) 
^ Jo 
where ji{t) = -xi{T)dxiV{xi+i{T) - xi{t)) 

We take average over all local heat current to find 
the global heat current, Ji{N) = Y^^Ji^ J^/{N - I). 
The steady state local temperature is given by k^Ti — 
{xi{t)'^)t. First we simulate the model for temperature 
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FIG. 2: Plot of Ji(iV) with iV for the ID DW model with 
k2 — ~l, k4 = 1 and for the FPU chain with k2 — k^ — 1 
(Length up to 1024) at high temperatures, Ti, — 1 and Tr — 2. 
The dotted straight lines correspond to the asymptotic length 
dependence of Ji(N) for the two models at this length. The 
inset shows local temperature profiles of the ID DW and FPU 
models for A^ = 1024. In all plots, 7l = 7r = 1. 



of the heat baths being large compared to Ti. For ex- 
ample, we choose temperature of the baths are, Tl = 1 
and Tr = 2 for /ca = -1, k^ = 1 and ti = 0.25. We find 
that the asymptotic length dependence of heat current 
for the ID DW model at high temperatures is similar 
to the FPU chain which is k ~ N^^^. We plot length 
dependence of Ji{N) for the DW model in Fig[21 and 
compare it to that of the FPU chain for similar length 
scales. We find Ji{N) falls relatively slower with increas- 
ing length for the FPU chain compared to the ID DW 
model for the simulated length here. The asymptotic be- 
haviour K ^ N^^^ in FPU chain was shown for a chain 
length of order of 10^ [7|. The local temperature T; of 
both the models has been included in the inset of Fig[2] 
for A^ — 1024. It also shows better local thermalization 
for the ID DW model than the FPU chain for shorter 
lengths. 

Now we reduce the mean temperature of the heat baths 
below Tl. The steady state heat current at low temper- 
ature in the DW model shows faster decay with length 
compared to high temperature. It has been claimed in 
Ref. J| on the basis of Green-Kubo formula that the low 
temperature thermal conductivity is independent of N 
for Nose- Hoover heat baths. It is worthy to mention that 
the uncertainty in the simulated thermal conductivity for 
the DW model in Fig.l of Ref. [J] is quite large. Therefore 
it is hard to conclude from their results that the low tem- 
perature heat conduction is normal in the ID DW model. 
We simulate the ID DW model here with Langevin heat 
baths at different mean temperatures (see Fig[3]). We 
find that it is difficult to simulate longer chain at lower 
temperatures keeping the numerical errors in simulated 
heat current relatively small. Also the length dependence 
of heat current varies non-monotonically with mean tem- 
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FIG. 3: Plot of Ji(iV) with iV for the ID DW model 
with ^2 = —1, ^4 ~ 1. The dotted straight Unes cor- 
respond to the asymptotic length dependence of Ji{N) 
at different mean temperatures. The inset shows a — 
Log[J(128)/J(256)]/Log[2] with T = (Tl + Tr)/2 where 
(Tr - Tl) < r. In all plots, 7l = 7r = 1. 



perature for a finite length as shown in the inset of FiglS] 
The inset shows the exponent a in Ji{N) ~ N~°' where 
N is varied between 128 and 256. Thus it becomes diffi- 
cult to find the asymptotic exponent of thermal conduc- 
tivity at low temperatures unambiguously. We are not 
able to get normal heat conduction in the ID DW model, 
but find that asymptotic length dependence of heat cur- 
rent varies between J{N) ~ 1/iV"-™ to J{N) - 1/N°-^'^ 
at different mean temperatures. We plot local temper- 
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FIG. 4: Plot of local temperature Ti with scaled site l/N for 
the ID DW model (fc2 = -1, fc4 = 1) and the FPU chain (fc2 = 
^4 = 1) for different chain length N at low temperatures. The 
inset shows local temperature profiles of the ID DW model 
for large A'^. In all plots, 7l = 7r = 1. 

ature Ti with site I of the ID DW model for different 
lengths, and compare them with that of the FPU chain 
at low temperature (see FigH]). While the ID DW model 
shows local equilibration even for smallest length con- 



sidered, there are large jumps in Ti at the boundaries of 
the FPU chain for smaller lengths, and the jumps become 
smaller with increasing length. It confirms that it easy to 
get local thermalization for the DW inter-atomic poten- 
tial compared to the FPU type potential. We show local 
temperature profiles of the DW model for longer lengths 
in the inset of Fig|4l They show a nice linear profile 
across the sample which signals normal heat conduction. 
Thus, though we could not find length independent k, it 
seems that heat conduction is normal for ID DW model 
at low temperature as it was concluded in Ref. [4] . 

Next we study the DW model in 2D and 3D. We 
consider a d-dimensional hyper-cubic lattice with the 
DW-type inter-atomic potential. Here we denote the 
lattice points by the vector n — {ni, n2, ..., n^j} with 
7ii — 1,2,..., TV, n,j=2,..d — l,2,...,W. The scalar dis- 
placement of an unit mass particle at the lattice site n is 
given by x^ . The Hamiltonian of the lattice is given by 



H.^jZri 



X!^Dw(a 



-e), 



(5) 



where e denotes the 2d nearest neighbors of any site. We 
couple all the particles at layers ni = 1 and ni = TV to 
Langevin heat baths, at temperatures Tl and Tr respec- 
tively. We again use fixed boundary conditions along the 
direction of transport, namely the 1 direction, and em- 
ploy periodic boundary conditions in the other (d-1) di- 
rections. We write Langevin equations of motion for the 
particles coupled to the baths, these are similar to EqsH] 
The corresponding white noise terms are related by fluc- 
tuation dissipation relations in EqlH We then simulate 
these equations using the velocity- Verlet scheme and de- 
termine the steady-state heat current and temperature 
profile. The time averaged local heat current along the 
1 direction is defined following Eq|4l We take average 
of local heat current over each layer to get average heat 
current through a layer with constant ni consisting of 
W' = W^"^ particles. In the steady state the average 
heat current through each layer is equal, and stationar- 
ity can be checked by testing how accurately this equality 
is achieved. Finally we take average of heat current over 
all layers to find the steady-state heat current Jd{N) in 
the d-dimcnsional lattice Uj, llj| . 

Recently the dependence of heat current Jd(N) on the 
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width W of the FPU lattice has been studied [10|, [ll| . It 
has been shown that for a fixed length TV, the heat current 
Jd{N) falls with increasing W but saturates quickly to 
the true higher dimensional heat current. The crossover 
width Wc at which value a crossover from ID to higher 
dimensional behaviour takes place, increases slowly with 
TV. We plot J2(TV) for the 2D DW model with lu = W/N 
for different values of TV in the inset of Fig [5] We find that 
the crossover from ID to 2D nature occurs at decreasing 
cj for increasing TV. It indicates that we can calculate 
heat current accurately in higher dimensional longer lat- 
tices by simulating relatively smaller W . We employ this 
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FIG. 5: K = NJ2(N)/AT with iV on (a) top: linear-log plot 
(b) bottom: double linear plot for the 2D DW model at high 
and low temperatures with k2 — —1, ki — 1. The inset in the 
bottom shows J2{N) with the normalized width lj = W/N 
for Tl = 1, Tr = 2. In all plots, 7l = 7r = 1. 



observation in our simulations for longer lattices in both 
2D and 3D. 

Next we check the dependence of heat current J2{N) 
on TV for the 2D DW lattice at high and low tempera- 
tures. An estimate of the crossover temperature T^ sep- 
arating two different structural phases in higher dimen- 
sions is more complicated than ID, but we can roughly 
estimate it as T^ ~ d Ti for a d-dimensional homogeneous 
lattice. The behaviour of heat current at high temper- 
ature of the baths, Tl = 1 and Tr = 2 is shown in the 
top of FiglS] using linear-log scale which clearly shows 
K ^ log(A^). Thus we confirm the predicted logarithmic 
behaviour of 2D FPU lattices in our simulations for the 
DW model with relatively shorter lengths. It is inter- 
esting to mention that this length dependence of k in 
2D DW model at high temperature is different from the 
previous simulated behaviour of 2D FPU lattices. The 
variation of the temperature gradient is nonmonotonic as 
a function of distance across the lattice in 2D at high tem- 
peratures (see Figl?]) which indicates non-diffusive trans- 
port as argued in Ref.pjJ. The length dependence of 
steady state heat current J2(iV) at low temperature is 
nonmonotonic with increasing mean temperature which 



is similar to the ID case. We find clear signature of nor- 
mal heat conduction in the 2D DW model at low tem- 
perature as shown in the bottom of FigO using double 
linear scale. The corresponding local temperature profile 
is fully linear. 
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FIG. 6: K = NJ3{N)/AT with N for the 3D DW model at 
high and low temperatures with ^2 = —1, ^4 = 1. In all plots, 

7l = 7r = 1- 

Finally we consider behaviour of J3 {N) for the 3D DW 
model. We readily find normal diffusive heat conduction 
at low temperature for quite shorter lengths. This is 
shown in Figl6] using double linear plot for bath tem- 
peratures Tl = 0.4 and Tr — 0.6. The corresponding 
local temperature profile is also linear as shown in the 
inset of Fig[71 The error bars in our simulated heat cur- 
rent for longer lengths are significant (about 2.5% for 
N=128), thus we could not confirm unambiguously nor- 
mal heat conduction for the 3D DW model at high tem- 
perature. Though we find signature of normal heat con- 
duction showing length independent k at high tempera- 
tures (shown in FigH] for Tl = 1.0 and Tr = 2.0). We 
strongly believe that heat conduction in 3D at high tem- 
perature for the DW model is normal and simulation for 
longer lengths would confirm it (even for much smaller 
length compared to Ref.[ll|)). We justify our claim for 
normal heat conduction for the 3D DW model by com- 
paring local temperature profiles of ID, 2D and 3D DW 
models at high temperature. It is shown in FiglT] which 
shows that the variation of the temperature gradients in 
both ID and 2D are nonmonotonic as a function of dis- 
tance across the lattice while for the 3D DW model it is 
almost monotonic and linear which are characteristics of 
normal diffusive behaviour 
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A dimensional crossover of thermal transport is re- 
cently observed in few- layer graphene [l5|. The type of 
crossover in thermal transport discussed here can be real- 
ized for the case of displacement ferroelectric transitions 
in which a crystal undergoes a structural transition from 
high temperature paraelectric state with higher symme- 
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FIG. 7: The plot of local temperature of ID, 2D, 3D DW 
model for TL = 1, Tr = 2. The inset shows local temperature 
for 3D DW model for Tl = 0.4, Tr = 0.6. In all plots, 
7l = 7r = 1- 



try to a low temperature ferroelectric state with a barrier 
at the normal lattice site |16| . 

The author is indebted to Prof. A. J. Sievers for valu- 
able discussions. 
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